Mon. Not. R. Astron. Soc. 000, [I]-?? (2006) Printed 5 February 2008 (MN IATrX style file v2.2) 



Is there a caustic crossing in the lensed quasar 
Q2237+0305 observational data record? 
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ABSTRACT 

We re-investigate the gravitationally lensed system Q2237+0305 data record to quan- 
■ tify the probability of having a caustic crossing in the A component. Several works 

\^ | assume that this is the case, but no quantitative analysis is available in the literature. 

. We combine the datasets from the OGLE and GLITP collaborations to accurately 

trace the prominent event in the lightcurve for the A component of the system. Then 
the observed event is compared with synthetic light curves derived from trajectories in 
magnification maps. These maps are generated using a ray-tracing technique. We take 



Of 



more than 10 9 trajectories and test a wide range of different physical properties of the 
, lensing galaxy and the source quasar (lens transverse velocity, microlens mass, source 

intensity profile and source size). We found that around 75% of our good trajectories 
(i.e. that are consistent with the observations) are caustic crossings. In addition, a 
high transverse velocity exceeding 300 km s _1 , a microlens mass of about 0.1 M© and 
a small standard accretion disk is the best parameter combination. The results justify 
r> \ the interpretation of the OGLE-GLITP event in Q2237+0305A as a caustic crossing. 

Moreover, the physical properties of the lens and source are in very good agreement 
with previous works. We also remark that a standard accretion disk is prefered to those 
simpler approaches, and that the former should be used in subsequent simulations. 
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1 INTRODUCTION 

The system QSO 2237+0305 is formed by a distant source 
quasar at zq — 1.695 and a spiral lens galaxy at zq = 0.039 
that quadruples the images of the source (Huchra et al. 
1985). Nowadays, it is one of the best studied gravitation- 
ally lensed quasar systems which has been monitored by a 
number of groups in the last years (Corrigan et al. 1991, 
0stensen et al. 1996, Vakulkik et al. 1997, Wozniak et al. 
2000, Alcalde et al. 2002, Schmidt et al. 2002). Among these 
monitoring campaigns, two of them have produced excelent 
light curves: OGLE (Wozniak et al. 2000) and GLITP (Al- 
calde et al. 2002). The former collaboration covers a time 
frame of several years, while the latter a few months but with 
a daily sampling, bad weather or technical problems aside. 
Interestingly, the two datasets overlap each other during a 
period of time. 



The recent datasets have been used to investigate sev- 
eral physical properties of the system. Limits on the size of 
the quasar emission region and the mass of the central black 
hole have been placed by Wyithe et al. (2000a), Yonehara 
(2001), Shalyapin et al. (2002), Goicoechea et al. (2003) and 
Kochanek (2004); the mass range of the microlenses in the 
lensing galaxy has been established by Wyithe et al. (2000b) 
and Kochanek (2004); and limits on the transverse velocity 
of the lensing galaxy have been reported by Wyithe et al. 
(1999) and Gil-Merino et al. (2005). To explain the promi- 
nent OGLE-GLITP event in the brightest component (A), 
some of these authors have assumed a caustic crossing (e.g., 
Goicoechea et al. 2003, Moreau et al. 2005). This interpreta- 
tion of the OGLE-GLITP event in the A component (caus- 
tic crossing) is supported by the best trajectories in terms 
of \ 2 that appear in Kochanek (2004). Although Kochanek 
(2004) also showed other relatively good trajectories passing 
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Figure 1. V-band photometry of the A component of 
Q2237+0305 for both data sets, OGLE (open circles) and GLITP 
(filled circles). There is still no calibration between the two sets. 

through complex magnification regions, the caustic cross- 
ing scenario became popular for studying the nature of the 
source quasar, e.g., reconstruction of the intrinsic bright- 
ness profile or estimation of the size ratios (Bogdanov & 
Cherepashchuk 2004, Goicoechea et al. 2004, Koptelova & 
Shimanovskaya 2005). 

In this contribution we merge the available data of the 
prominent OGLE-GLITP event (it was monitorized by both 
OGLE and GLITP collaborations) , building a single dataset 
for the A component. We try to find the probability of hav- 
ing a caustic crossing in this component considering different 
source profiles and sizes, microlenses masses and lens trans- 
verse velocities. In this way we both justify those works that 
assumed approximations to caustic crossings and future as- 
sumptions. 

The paper is organised as follows: Sec. [5] describes the 
way in which the datasets from the different OGLE and 
GLITP collaborations are merged to produce a single set; 
Sec. |21 explains the method used to produce the synthetic 
light curves; Sec 2] establishes the criteria to compare the 
synthetic light curves and the observational record, and gives 
the results of these comparisons; finally, in Sec.|S]we include 
our conclusions. 



2 DATA MERGING 

We are interested in the A component of Q2237+0305. The 
system was monitored by the phase II of the OGLE 1 project. 
During part of that period, the GLITP 2 project obtained 
optical frames of the same lens system. In particular, the 
overlapping period was approximately between the Julian 
days 2451400-2451600, corresponding to from October 1999 
until February 2000. The corresponding observational data 
set from both campaings in the V filter is displayed in Fig. El 
From now on we use a simpler scale of time: JD-2450000, 
so a peak of flux around day 1500 is seen in Fig. 1. 

In order to merge the OGLE and GLITP data sets, the 

1 http: / / ast rouw.edu.pl/ ~ogle/ 

2 data from [http: / /www.iac.es/proyect/gravs-le ns / GLITP / 
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Figure 2. OGLE (open circles) and GLITP (filled circles) data 
sets have been merged following a "re-reduction" of the original 
OGLE data and expressions in Equation [3] (see text for further 
details). 

first problem we found was the night sampling. During the 
OGLE campaign, two data points were consecutively ob- 
tained every night. The initial step is thus to "re-reduce" 
these data points to get one per observational day. In this 
way, when two points are available for a given night, the 
observational date and flux will be the mean of the indi- 
vidual observational dates and fluxes. The estimation of the 
errors for the new points has two cases: if the individual 
error bars overlap, the error is the mean of the individual 
errors; if they do not overlap, then the error is estimated by 
the square root of the sum of the square of the average of the 
individual errors plus the square of the difference between 
the mean flux and one of the two individual fluxes. 

Once we have the same night sampling for both data 
sets, the next step is to merge them. To do that, we try to 
find a linear expression that relates the data sets: 

F GL E (t)=a + b- FGLITp(t) (1) 

where F(t) are the fluxes and a and b are the constants to 
be determined. We treat the problem as a x 2 minimization 
(e.g., Ofek & Maoz 2003, Bogdanov & Cherepashchuk 2004). 
The best fit found is 

F OG LE(t) = -0.097 + 1.183 • F GL iTp(t) (2) 

In the GLITP photometric system, the "new" OGLE 
data points are then 

F'ogle (t) = 0.082 + 0.845 • F OG L E (t) 

v'ogle = 0.845 • oogle (3) 

The resulting lightcurve is presented in Fig. [5] This is 
the dataset we will compare to simulations. 



3 SIMULATIONS 

Prior to produce synthetic light curves, we first build magni- 
fication maps where these will be drawn. The 2-dimensional 
magnification maps are done using a ray tracing method: 
light rays are back-traced from the observer to the source 



plane through a distribution of mass sitting on the lens plane 
(Kayser et al. 1986, Wambsganss 1990). The distribution of 
mass responds to a macromodel of the system, that gives the 
two relevant parameters: the convergence k or the distribu- 
tion of mass in compact objects at the position of the beam 
and the shear 7 or the influence of matter outside the beam. 
In general, lightcurves of lensed quasars might contain both 
microlensing and instrinsic variability. So in order to analyse 
microlensing alone one should correct for the time delay and 
then subtract to components. However, in the case of Q2237 
the time delay is expected to be of a few hours (Vakulik et 
al. 2006). In addition, in the case of the period analysed 
here for Q2237, the flatness of component D indicates that 
intrinsic variability does not contribute significantly to the 
behavior of the components (see Shalyapin et al. 2002). We 
use the parameters for the A component given by Schmidt et 
al (1998): k = 0.36 and 7 = 0.40 (we note that in Gil-Merino 
et al. 2005 slightly different models where indistinguishable, 
driving to the same conclusions). The resulting density of 
rays at a point in the source plane is proportional to the 
microlensing magnification of the source at that position. 
The space scale factor is the Einstein radius, defined in the 
source plane as 
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TE = 
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where M M i ons is the mass of the microlenses and Di are 
angular distances: D s is observer-source, Dd observer-lens 
and Dds lens-source. Setting the space scale (te), the value 
of the microlens mass is directly related to the choice of 
the cosmological model. The magnification maps are then 
convolved with a particular intensity profile for the source, 
which is characterized by a typical radius. We must also as- 
sume a certain effective transverse velocity for the source 
(transverse displacement of the source in the source plane 
divided by elapsed time in the observer's clock) Vt = — (1 + 
zg) _1 (Ds I D d )v t> i ens (see e.g. Gil-Merino et al. 2005 and 
references therein for the relationship between the effective 
transverse motions of the source and the lens). The syn- 
thetic microlensing light curves are obtained extracting lin- 
ear trajectories across the magnification patterns, since once 
all physical parameters are fixed and knowing the observa- 
tional baseline of 200 days, we can calculate the length of 
these trajectrories. We note here that we do not take into 
account the additional effect of the proper motion of the mi- 
crolenses. This effect would play a significant role in large 
time-scale lightcurves but in short time-scale campaigns the 
magnification maps can be considered static to this respect. 
Including the proper motion of microlenses, we would in- 
crease the number of caustic crossings in any case, so the 
results here are conservative. 

To carry out the simulations we take H = 66 km s _1 
Mpc -1 . For a given value of te in Eq. (4), the microlens 
mass depends on the matter-energy density of the universe. 
We focus on both the Einstein-de Sitter model (EdSM: Qm 
— 1, Qa = 0) and the current concordance model (CM: 
Q.M — 0.3, S1a = 0.7), so each cosmology leads to a differ- 
ent microlens mass, but the final statistics on the caustic 
crossings remains the same. Taking a particular value of Vt, 
we also have two different values of the effective transverse 
velocity for the lens, one corresponding to the EdSM and 
another associated with the CM. In Tablc^all the physical 



intensity 


source 


Aliens (Mq) 


v t ,i ens ( km s *) 


profile 


radius (cm) 


EdSM/CM 


EdSM/CM 


Uniform 


6-10 14 


0.05/0.04 


100/75 


Gaussian 


2-10 15 


0.10/0.08 


300/220 


Standard 


610 15 


0.6/0.5 


600/445 



Table 1. We tested three different intensity profiles, source radii, 
Einstein radii and effective transverse velocities for the source. 
The resulting simulations were obtained for all the possible com- 
binations between these parameters. Here, standard profile means 
intensity profile of a Newtonian geometrically-thin and optically- 
thick accretion disk. The source radius is the typical radius of the 
intensity profile, so it is the radius of the quasar for a uniform 
disk, but is smaller than the quasar radial size (containing 95% 
of the total brightness) in the other cases (see Shalyapin et al. 
2002). On the other hand, the values of M^i^g and v t: i ens de- 
pend on the cosmological model, and we present their values for 
the EdSM and CM with H = 66 km s" 1 Mpc -1 (see main text). 



values are presented. We use a set of axisymmetric sources, 
including brightness distributions enhanced at the centre of 
the source (standard accretion disk and and Gaussian pro- 
files) as well as a uniform brightness distribution. Their typ- 
ical radii are consistent with previous studies by Wyithe et 
al. (2000a), Yonehara (2001), Shalyapin et al. (2002) and 
Kochanek (2004). Wyithe et al. (2000b) obtained that the 
most likely value for the microlens mass is in the range 
0.01-1 Mq, while Kochanek (2004) used a CM and reported 
an interval 0.003-0.1 M . The values of M M ; ens in Table 1 
(EdSM and CM) vary between a lower limit slightly above 
0.01 Mq and an upper limit slightly below 1 Mq, and they 
reasonably agree with previous work. With respect to the 
lens transverse velocities, Wyithe et al. (1999) claimed that 
vt,iens < 500 km s _1 (EdSM). Moreover, using an interme- 
diate microlens mass of 0.1 M , Gil-Merino et al. (2005) 
also inferred a constraint v t ,iens < 630 km s _1 (CM). Our 
values of Vt,i en s in Table 1 are basically consistent with these 
dynamical constraints. 

Taking into account that each convolution gives a differ- 
ent map, we have produced a total number of 108 magnifica- 
tion maps, 2x2 re on 2000 pixels a side each. Considering all 
possible combinations of the transverse velocities, the num- 
ber of time maps (by transforming the spatial axes into time 
axes) is 3 times greater. The number of trajectories tested 
on each magnification pattern is 10 7 . The total CPU time 
for these simulations, excluding the generation of maps, was 
around 720 hours on a Pentium 4. 

But could Nature be conspiring against us? It could 
happen that the variability detected in the A image (see 
Fig. |2J is in fact a cusp passing and what we see is a com- 
bination of a symmetric microlensing event due to the cusp 
plus a gradient originated by intrinsic variability. If this is 
true, that gradient should also appear in component D dur- 
ing the same period, because the time delays are of only 
the order of hours. Since this component is flat, the only 
explanation is that microlensing and intrinsic variability are 
mutually canceled in image D. Fortunately we can test this; 
we construct a magnification map for the D component of 
4500 pixels on a side covering a physical length of 10 Ein- 
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Af M ; ens (Mg) 


trajectories with 5^1 


EdSM/CM 




0.05/0.04 


26% 


0.10/0.08 


56% 


0.6/0.5 


18% 



Table 2. More than a half of the good trajectories suggests a mass 
of the microlenses of about 0.1 Mq, among the values tested. 

stein radii, and we then analysed the gradient of 10 6 tracks 
of a length equal to the 200 days period of observations, for 
each combination of source size and transverse velocity. We 
search for gradients between a 15% and a 23%, the ranges 
associated with the gradient seen in the A component. We 
found that for the biggest source size considered in this study 
only 0.15% of the cases showed a gradient in that range. The 
percentage was slightly bigger, 1.16%, for the smallest source 
size. These results were obtained for transverse velocities of 
600 km s _1 . Considering velocities of 300 km s -1 the per- 
centage was 0.5% for the smallest source size and 0.001% for 
the biggest one. Thus, we conclude that the possibility of a 
conspirancy between microlensing and intrinsic variability is 
extremely unlikely in this work. 



4 ANALYSIS AND RESULTS 

Each simulated lightcurve is compared to the observational 
one. The way we did it is a x 2 test. We estimate the differ- 
ence between each \ 2 value and the the number of degrees 
of freedom (u), and then we select the good fits verifying 

S = J^=^ < 1 (5) 

where 8 is the relative deviation of \ 2 (Shalyapin et al. 
2002). All the synthetic trajectories that satisfy Eq. arc- 
considered as physical processes consistent with the obser- 
vational lightcurve, i.e., good trajectories. Then we compute 
how many of them are in fact caustic crossings in their re- 
spective magnification maps. After generating a total num- 
ber of 3.34- 10 9 synthetic lightcurves, we found that 54 of 
them satisfied Eq. |5] with a (x 2 ) = l-125. Among them, 14 
were not caustic crossings, but cusp crossings. This means 
that the event seen in the A component of Q2237+0305 can 
be explained as a caustic crossing in a 74% of the good tra- 
jectories generated in this contribution. The three panels of 
Figure 3 show some good trajectories in our magnification 
maps. Caustic crossings appear in the top and bottom panels 
of Fig. 3. However, the middle panel includes less probable 
(1:4) cusp crossings. 

Interestingly, the results can also be sorted attending 
to different criteria. For example, all the trajectories match- 
ing Eq. correspond to a standard accretion disk. None of 
the other intensity profiles seem to be viable to explain the 
prominent event. Due to the size of the maps, we cannot 
properly test the extreme case of a standard source with the 
larger radius. However, considering the other two radii, the 
smaller one is clearly favoured (91% of the trajectories with 

5 ^ 1). In addition, more than a half of the trajectories with 




Figure 3. Caustic crossings (top and bottom panels) and cusp 
crossings (middle panel) that agree with the OGLE-GLITP 
event in Q2237+0305A (see Fig. 2). Their associated synthetic 
lightcurves verify Eq. [5] 

5^1 correspond to a microlens mass of about 0.1 Mq (Ta- 
ble and 87% of them to transverse velocities of 450-600 
km s _1 (Table 0J . 



5 CONCLUSIONS 

Several previous works using either OGLE or GLITP (or 
both) data sets assumed a caustic crossing occurring in the A 
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V t ,lens (km S 1 ) 


trajectories with S sC 1 


EdSM/CM 




100/75 


0% 


300/220 


13% 


600/445 


87% 



Table 3. The maximum transverse velocity tested is favored by 
the simulations against lower values. This is interesting because 
the previous strongest upper limits for the effective transverse 
velocity were of ~ 600 km s _1 . 

component of Q2237+0305 (see, e.g., Shalyapin et al. 2002). 
In this contribution we tested, via ray-tracing simulations, 
whether the event seen in component A running from Oc- 
tober 1999 to February 2000 is compatible with a caustic 
crossing and to which extend. From a very large number of 
simulated trajectories (> 10 9 ) we found that only several 
tens were compatible with observations. These good trajec- 
tories are related to two different kinds of physical scenarios: 
caustic crossing (see the top and bottom panels of Fig. 3) 
and cusp crossing (see the middle panel of Fig. 3) . A caustic 
crossing scenario has a probability of about 75% (40 good 
trajectories), whereas a cusp crossing scenario (do not con- 
fuse with a source passing close, but outwardly, to a cusp) 
has a smaller probability of about 25% (14 good trajecto- 
ries). We checked different lens transverse velocities with an 
upper limit of 600 km s" 1 (Wyithe et al. 1999, Gil-Merino 
et al. 2005). Surprisingly, small or intermediate velocities 
(^ 300 km s _1 ) do not seem to be the most probable ones, 
which means that the true transverse motion might be close 
to that upper limit. Our results also suggested that the mass 
of the microlenses at the location of image A should be very 
close to 0.1 Mq, which is in the centre of the range derived 
by Wyithe et al. (2000b) and is the upper limit reported by 
Kochanek (2004). Yonehara (2001) and Gil-Merino et al. 
(2005) also assumed this mass to obtain their strongest con- 
straint on the quasar size and the lens transverse velocity, 
respectively. Interesting is also the fact that the best results 
were obtained using a standard accretion disk as intensity 
profile and that the smallest source size tested was the fa- 
vored one. This last result is also in agreement with some 
previous work (Shalyapin et al. 2002, Kochanek et al. 2004). 
We note that we did not investigated different inclinations 
of the disk nor more sophisticated accretion disk models, 
which is an interesting issue considered as a second order 
effect but to be addressed in a future work. Very recently, 
Mortonson et al. (2005) found that the generic microlensing 
fluctuations are relatively insensitive to all properties of the 
source models except the half-light radius of the disk. Here 
and in other related papers it is showed the feasibility of a 
discrimination between different intensity profiles through 
the observation and analysis of caustic/cusp crossings. 
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